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BACKGROUND OF THE INVENTION 
[0001] The present invention relates to a method of keying a borehole in a multi-attributes 
seismic block, that it is to say it aims to obtain the best position of a drilling trajectory in the 
seismic block. 

State of the art 

[0002] The method according to the invention applies to seismic blocks. A seismic block is 
formed by the juxtaposition of seismic traces in a volume. The term "seismic section" refers to a 
slice of seismic block. A seismic section offers a view of the juxtaposition of the seismic traces 
contained in the cutting plane. In a seismic image, the luminous intensity of a pixel is 
proportional to the represented seismic attribute. 

[0003] The keying of a borehole in a seismic domain has already been the subject of some 
works. 

[0004] In particular, patent FR-A-2 737 909 proposes a method of positioning a drilling 
trajectory on a seismic section by the modelling in a synthetic sum domain of locations of the 
trajectory for which the pair (dip, azimuth) is known then by migration in time of the synthetic 
image. 

[0005] This method allows the bore to be placed on a seismic section by a deterministic 
calculation method. 

[0006] Other calculation methods have been proposed. However, they are all deterministic. 
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Contribution of the invention 

[0007] The approach of the method according to the invention is non-deterministic in nature. It 
is accepted that any measurement is marred by error and the intention is to determine the best 
possible keying by a scanning of all the keying possibilities in order to select the best of them. 

[0008] The method according to the invention determines the value of the keying through the 
convergence quality in learning phase of a set of neural networks. 

[0009] The novelty lies in the fact that the neural networks are not only used as universal 
estimators but also as measuring instruments. 

[0010] Thus each network produces a "distance" keying measurement and the best keying is 
determined by the minimum distance. 

[0011] This operation is analogous to the selection of the minimum distance from a set of 
distance measurements taken with any measuring instrument, such as a graduated rule. 

Definition of the invention 

[0012] The subject of the present invention is a method of keying, in a space presenting two 
spatial dimensions and one temporal dimension, a signal S measured in positions U subject to an 
uncertainty, on a set of N signals measured in determined positions, the N +1 signals having their 
temporal origin in a same plane, the said method involving: 

- re-sampling the N + 1 signals in order to place them all in an identical 
sampling range, 

filtering the signal S in order to place it in a range of frequencies that is identical 
to that of the N signals, 
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and characterized in that it also involves: 

defining for each position U associated with the measurements of the signal S a 
same neighbourhood of places V in the spatia-temporal space centred on the position U, 

producing a layered neural network RN V for each location V in the neighbourhood 
of U, each network having an entry vector of dimension N associated with the measurements of 
the N signals and a scalar exit associated with a measurement of the signal S, 

- for each neural network RN V , defining a learning set such that the entries are the 
collection of all the vectors of measurements of the N signals situated at locations V and the exits 
are the collection of the values of the signal S at positions U for all the positions U, 

fixing a predetermined number of iterations Nit for all the neural networks and 
launching the learning phases of all the networks, 

for each neural network RN V , calculating the value of the integral £ v of the 
function giving the error committed by the network at each iteration, from iteration 1 to iteration 
Nit, 

for each surface spatial position V k of the neighbourhood with coordinates 
(Xk,Yk,tO). selecting in the time dimension the pair of locations Vl k (x k ,y k ,to,), V2 k (x k ,y k ,t 2 ) of 
the neighbourhood which correspond to the two smallest local minima of the two integrals (£ vl k , 

r 2 k ), 

- for each surface spatial position Vk of the neighbourhood, retaining from among 
the two positions VI k (x k ,y k ,to,) and V2 k (x k ,y k ,t 2 ) the position V m , for which the signal estimated 
by the respective neural networks RN vl k and RN v2 k presents a maximum variance, choosing from 
among the positions V m the position V ca] for which the integral £ v m is minimum. 
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Figures 

[0013] Figure 1 represents an electric signal positioned in a seismic block. 

[0014] Figure 2 illustrates the result of the filtering of an electric signal which adjusts its range 
of frequencies to that of the seismic signals. 

[0015] Figure 3 represents a process diagram showing the functioning of the neural networks in 
the method according to the invention. 

[0016] Figure 4 shows four error functions associated with the learning phases of four layered 
networks. 

[0017] Figure 5 illustrates the principle of selection of the best time drift. 

[0018] Figure 6 presents a neural network in which the entries have come from a time window 
taken on the seismic signals. 

[0019] Figure 7 shows the result of the estimation of an electric recording realized by the neural 
network and applied to the best keying point. 

Complete description 

[0020] The method according to the invention is a method for keying a borehole in a seismic 
block. With reference to the figures, an embodiment of this method is given below. 

[0021] A seismic acquisition permits the production of a seismic block the purpose of which is 
to represent the volume of subsoil through which the seismic wave has travelled. A seismic block 
is constituted by a set of seismic traces that have been summed, that is to say pre-processed in 
order to obtain the best possible view of the underground reflectors. The reflectors are interfaces 
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between two media presenting different speeds for the seismic wave. The reflectors generally 
represent limits of geological layers. 

[0022] Let B be a seismic block. The seismic block is defined by a three-dimensional matrix M. 
Two dimensions of the block B are linked to two horizontal spatial axes, commonly called (Ox, 
Oy). The third dimension of the block B is linked to a vertical temporal axis, commonly called 
OL. 

[0023] As the seismic block is an image of the subsoil, each element of the matrix M is 
associated with a point with coordinates (x b i, yw, z b i) in the subsoil. The time component t b i is 
linked to the depth zbj by a relationship of the type: 

zbi = v(tbi) ' tbi 

[0024] where v(tbj) is a function which gives the speed of the seismic wave as a function of the 
depth. 

[0025] Each element of the matrix M contains a seismic attributes vector (An, Ai 2 . . . A in ) at the 
point with coordinates (xbj,ybi, tbj). For example, the attributes (An, A i2 ...Ai n ) measure the 
energy of the seismic wave, also called seismic amplitude. 

[0026] The term offset refers to the distance which separates the transmitter-receiver pair 
associated with a seismic trace. A seismic amplitude measurement is then determined by the 
summation of the seismic amplitude measurements according to predetermined classes of offsets. 
Typically, the traces will be summed according to a division into three classes of offsets: near 
offsets, far offsets and very far offsets. 
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[0027] Let F be a borehole drilled in the volume of subsoil covered by seismic block B. The bore 
F travels through the subsoil along a drilling trajectory defined by a set of points with spatial 
coordinates (xfj, yfj, zfj). 

[0028] Any oil exploration drilling gives rise to the recording of numerous electric 
measurements in the depth of the bore. For each point with coordinates (xfj, yfj, zfj) of the 
trajectory of the borehole, a set of electric measurements (Ei, E 2 , . E n ) is thus acquired. 

[0029] The dimension of the depths zfj, is the unit of length; the measurements are taken in 
metres or in feet. A first operation to place the bore in the seismic block is to transform the 
length value into a time value. This transformation is easily carried out with the help of 
relationship (1) since the rates of propagation of the seismic wave in the subsoil in the 
neighbourhood of the bore can be measured, these rates being measured directly in the borehole. 

[0030] The result of the conversion is a temporal bore trajectory defined by a set of points U with 
spatial coordinates (xfj, yfj, zfj). 

[0031] The bore is then locked into the seismic block by the initial coordinates of its trajectory or 
surface coordinates (xf 0 , yf 0 , tf 0 ). 

[0032] Figure 1 shows schematically a portion of a seismic block (seismic cube) presenting steps 
in space dx and dy in the two horizontal dimensions and a step in time dt in the vertical 
dimension. An electric signal (Log) is placed in the portion of the block. 

[0033] As electric recordings in boreholes are generally carried out every foot and the resolution 
of a seismic signal, a function of its pass-band, is of the order of a few metres, it is advisable in a 
first step to re-sample the electric and seismic signals in such a way that their sampling steps in 



NY02:518494.1 



11 



A36002- 072853.01 10 
PATENT 

the time dimension are similar. To do this, all of the electric measurements of the bore Fare re- 
sampled and all of the seismic amplitude measurements of the block Bare oversampled. In a 
second step, a filter is applied to the electric signals in order to bring them into a range of 
frequencies that is identical to that of the seismic signals. 

[0034] Figure 2 shows an example of filtering of an electric measurement. The left-hand 
rectangle called "Log" represents the amplitude of an electric signal sampled in time after 
application of the transformation (1). The right-hand rectangle called "Filtered Log" represents 
the electric signal the frequencies of which have been brought into the domain of the seismic 
frequencies by application of high-pass and low-pass filters. The top rectangle called "Frequency 
Spectrum", the horizontal axis of which graduates frequencies and the vertical axis amplitudes, 
shows the two frequency spectra associated with a seismic trace (Near) and with the filtered 
electric signal (Filtered Log). It can be seen that the electric signal which originally contained 
many more high frequencies than the seismic trace is now brought into the frequency domain of 
the seismic signal. 

[0035] There is thus defined for each point U constituting the trajectory of the bore and 
associated with an electric measurement Ej a spatia-temporal neighbourhood in the seismic block 
B in the form of a parallelepiped centred on the point U. Taking the seismic block B as an 
example, the parallelepiped has two horizontal spatial sides and one vertical time dimension, 
orthogonal to the spatial side. Let Dvx be the number of locations of the neighbourhood on the 
axis Ox, Dvy the number of locations of the neighbourhood on the axis Oy and Dvt the number 
of locations of the neighbourhood on the axis Ot. 
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[0036] There will be chosen, for an example, a dimension of the parallelepiped of 5 samples 
either side of the point U in each dimension. In this case Dvx, Dvy and Dvt all take 1 1 as a value 
and the neighbourhood is defined as a cube containing 1331 (1 1 x 1 1 x 1 1) seismic measurement 
points. The same neighbourhood applies to all the points of the trajectory, that is to say by 
extension to the trajectory itself. The neighbourhood is defined by a set of points V, The number 
of points included in the neighbourhood is called Npv. Npv is equal to the product Dvx x Dvy x 
Dvt and equals 1331 in the example above. 

[0037] Npv layered neural networks are then produced. For a cube with 1331 points, 1331 neural 
networks are thus produced, The entry vector of each network is a vector of seismic amplitude 
measurements each component of which belongs to a different class of offsets. The exit of the 
network is reduced to a scalar representing an electric measurement Ej. 

[0038] Figure 3 provides a process diagram, In this example, the neural networks have two 
entries and one exit. The left-hand rectangle called "Near" shows the values taken by the first 
entry of the neural network in the learning set. The central rectangle called "Far" shows the 
values taken by the second entry of the neural network in the learning set. The right-hand 
rectangle called "Log" shows the values taken by the exit of the neural network in the learning 
set. An example is thus constituted by two entry values taken at a determined time (Time) on the 
"Near" and "Far" signals and an exit value taken at the same time on the _'Log" signal. The 
ellipse called "RN" is a symbolic representation of the neural network between entries and exit. 

[0039] The purpose of the neural networks is thus to discover a relationship between the electric 
measurement EjU and the seismic measurements (A/, A 2 V , ... A n v ), that is to say a function f 
such that: 
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E, u = f(Ar,A 2 v ,...A n v ) 

[0040] A measurement superscripted by u indicates that its value is taken on the trajectory of the 
bore and a measure superscripted by v indicates that its value is taken in the neighbourhood. 

[0041] The neural networks all present the same architecture, that is to say the same number of 
hidden layers, the same number of neurons per layer, the same number of connections and the 
same initial conditions. 

[0042] A position V in the neighbourhood is chosen. The neural network RN V associated with 
the position V is used as follows. For all the points U of the trajectory of the bore, a learning set 
is produced for the network RN V by collecting all the vectors of seismic amplitude measurements 
situated at position V. The learning set for the neural network RN V is thus constituted by as many 
examples as there are points U on the trajectory. 

[0043] A predetermined number of iterations Nit is chosen and the learning phase of the neural 
network is launched by retro-propagation of the gradient. 

[0044] When the Nit iterations are effected, the value of the integral £ v of the function giving 
the quadratic error committed by the network on each iteration, from iteration 1 to iteration Nit is 
calculated. 

[0045] A new position V in the neighbourhood of the trajectory is then chosen and the previous 
steps are carried out with the neural network RN V| associated with the position V in order to 
obtain a new integral £ v . 



NY02:5 18494.1 



14 



A36002- 072853.01 10 
PATENT 

[0046] The operation is thus repeated until all the points of the neighbourhood have been 
scanned. For a cube with 1331 points, 1331 integral values are thus obtained. 

[0047] With a machine, these operations can be carried out sequentially with the same neural 
network reinitialized with each change of point in the neighbourhood. However, it is 
recommended to use this technique in parallel in machines with multi processors organized in 
calculation collectives, producing as many neural networks as the capacity of the machine 
configurations allows. The ideal would be to be able to obtain as many neural networks as there 
are points in the neighbourhood. 

[0048] All the points of the neighbourhood having been travelled through, for each surface 
spatial position Vk of the neighbourhood with coordinates (x k ,y k , t 0 ) to the number of Dvx x 
Dvy, there is retained in the time dimension the pair of locations Vl k (x k , y k , ti), (V2 k (x k , y k , t 2 )) 
of the neighbourhood which correspond to the two smallest local minima of the two integrals 

G vi k ,r 2 k ). 

[0049] Figure 4 shows four error functions associated with the learning phases of four identical 
neural networks. The horizontal axes are graduated in iterations and the vertical axes give the 
quadratic errors committed by the networks during their estimation of the electric signal. The 
values of the four integrals associated with the networks are equal to the hatched surfaces. In this 
example, these are the networks identified by the triplets (dx k , dy k l, dti) and (dx k , dy k , dt 2 ) which 
present the two smallest local minima of the integrals. 

[0049] The variance of the signal estimated by the neural networks RN vl k and RN v2 k is then 
calculated and the position V m for which the variance is maximum is retained in the pair 
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VI k (xk, yk, tj), (V2 k (x k , y k , t 2 )). Thus a set composed of Dvx x Dvy locations V m in the 
neighbourhood is defined. 

[0050] Figure 5 shows how the selection of the position V m is made from the positions V k . The 
upper rectangle gives the values of integrals as a function of the shift on the temporal axis for a 
chosen V k . In this example, the temporal neighbourhood is fixed at 20 temporal samples either 
side of a point of the trajectory. The drifts DEI and DE2 correspond to the two smallest local 
minima of the integrals. The lower rectangle represents the values of the variances of the two 
signals estimated by the neural networks associated with the drifts DEI and DE2. The drift DEI 
is retained as the best keying V m for the surface position V k . 

[0051] Finally, from among all the Dvx x Dvy locations V m , the position V ca i for which £ v m is 
minimum is chosen. The position V ca i is the position of the neighbourhood which ensures the 
best keying of the electric signal on the seismic signals. 

[0052] Another embodiment involves not informing the neural networks by a single 
measurement per class of offsets, but by a set of M values taken in a time sampling window on 
the trace, which window is centred on a point V of the neighbourhood. A time window equal to 8 
samples either side of the point V, that is to say 17 time values, will be taken for example. 

[0053] Figure 6 presents such a network for which the chosen time window is 2 samples either 
side of the point V for reasons of simplicity. The network thus has 5 entries per class of offsets. 
[0054] The figure shows two classes of near (Near) and far (Far) offsets which generate 2x5 
entries. 
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[0055] This modifies the architecture of the network which thus comprises an entry vector the 
dimension of which increases from NtoNxM. The exit of the network remains scalar. 

[0056] The method according to the invention then applies in identical fashion and leads to the 
selection of the best point of the neighbourhood as the best keying possible. This windowing 
technique has the advantage of permitting the network to carry out a deconvolution on the basis 
of the duration of the seismic pulse transmitted in the subsoil. The quality of the results is 
thereby greatly improved. 

[0057] Figure 7 shows the result of the estimation of the electric recording (Prediction) realized 
by the neural network from the seismic signals and applied to the best keying point. The dotted 
line (bore) represents the electric signal that is effectively measured. The very great quality of the 
estimation, which affirms the quality of the keying determined by the method according to the 
invention, will be noted. 
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